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Superfluidity and superconductivity have been studied widely since the last century in many 
different contexts ranging from nuclear matter to atomic quantum gases. The rigidity of these 
systems with respect to external perturbations results in frictionless motion for superfluids and 
resistance-free electric current in superconductors. This peculiar behaviour is lost when external 
perturbations overcome a critical threshold, i.e. above a critical magnetic field or a critical current 
for superconductors. In superfluids, such as liquid helium or ultracold gases, the corresponding 
quantities are critical rotation rate and critical velocity, respectively. Enhancing the critical values 
is of great fundamental and practical value. Here we demonstrate that superfluidity can be achieved 
for flow above the critical velocity through quantum interference induced resonances. This has far 
reaching consequences for the fundamental understanding of superfluidity and superconductivity 
and opens up new application possibilities in quantum metrology, e.g. in rotation sensing. 


The breakdown of superfluidity and superconductivity 
above the critical velocity and the critical current [1, 2], 
respectively, is caused by the production and growth of 
excitations. In a superfluid, when its flow velocity, or 
equivalently the velocity of a defect dragged through the 
fluid, exceeds the critical velocity Vc the creation of exci¬ 
tations becomes energetically favourable. This destroys 
the frictionless motion, as shown in experiments with su¬ 
perfluid helium [3] , and ultracold bosons [4] and fermions 
[5] . The critical velocity is defined as the maximum veloc¬ 
ity below which there is no (or more precisely a bounded) 
production of excitations. While (essentially) no excita¬ 
tions are present for subcritical velocities, a fast onset 
of growing excitations occurs for supercritical velocities 
that gradually decreases for further increased velocity un¬ 
til the kinetic energy of the fluid becomes so high that it 
dominates all other energy scales and any defects become 
unimportant again (Fig. f). 

In general, the rate of excitation growth depends on 
the microscopic details of the superfluid and on how it 
is coupled to the environment. The details of this cou¬ 
pling determine the dissipation mechanism causing the 
creation of excitations and ultimately the critical veloc¬ 
ity Vc- A simple way to estimate Vc and give a qualitative 
explanation of the breakdown of superfluidity is provided 
by the Landau criterion [6]. 

The Landau criterion is a cornerstone of our under¬ 
standing of the dynamical behaviour of superfluids, stat¬ 
ing that a superfluid flow is sustained against external 
perturbations or defects up to a critical value of the ve¬ 
locity [6]. Its elegance, power and usefulness rely both 
on simplicity and generality: there is no need to know 
the specific nature of the perturbation or the character¬ 
istics of the defects, no need to know the microscopic 
details of the superfluid, and no need to compute the ex- 



FIG. 1. Excitations in a superfluid as a function of veloc¬ 
ity of a defect dragged through the fluid at a given time. 
For velocities below a critical velocity (v < Vc), the produc¬ 
tion of excitations is suppressed and superflow persists: at 
longer times the amount of excitations remains very small. At 
higher velocities (v > Vc) a sharp onset of excitations (grow¬ 
ing with time) destroys the superfluid properties and only at 
very high velocities {v S> Vc) the kinetic energy becomes so 
high that the defect hardly affects the flow. This standard 
picture (dashed line) has to be adjusted when resonant quan¬ 
tum interference reinstates superfluidity with fully suppressed 
excitation growth at a series of discrete supercritical veloci¬ 
ties. The solid line shows an example for a rectangular defect 
shape (for the same plot at different times see Fig. 2). 

citation spectrum of the moving system; only knowledge 
of the low-energy excitation spectrum e(p) of the system 
at rest is required. Briefly, by applying a Galileo trans¬ 
formation to the co-moving frame it can be shown that 
for V < vl (where vl = minis the Landau criti¬ 
cal velocity [6]) the production of elementary excitations 
is energetically unfavoured. From the Landau criterion 
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it follows that for short-range weakly interacting Bose 
gases, the critical velocity is equal to the sound velocity 
c, so that here “supercritical” means “supersonic”. Su¬ 
perfluidity in a weakly interacting Bose gas flowing at 
a velocity greater than the Landau critical velocity was 
studied in [7]. Although the detailed analysis of differ¬ 
ent superfluid systems, including helium and ultracold 
gases, shows that the Landau criterion often quantita¬ 
tively overestimates the actual critical velocity Vc, espe¬ 
cially in the two- and three-dimensional case (see also 
a recent paper on one dimension [8]), the identification 
of a critical velocity above which the production of ex¬ 
citations destroys the superfluid motion is a criterion of 
paramount clarity and relevance. 

The Landau criterion is based on a perturbative treat¬ 
ment of “small” defects affecting the superfluid motion. 
The possibility to explore superfluid motion in a non- 
perturbative regime of parameters has attracted con¬ 
siderable interest, e.g. in non-perturbative and/or exact 
studies of the dynamical propagation of a superfluid in 
presence of “non-small” defects (of tunable shape and in¬ 
tensity) or periodic potentials [9]. The point we address 
in this paper is the very surprising possibility of stable su- 
percritical propagation (i.e. above the critical velocity) at 
a non-perturbative level. In the following we show that 
it is possible to have a set of velocities larger than the 
critical velocity Vc where superfluid flow exists and the 
production of excitations is completely suppressed. The 
key mechanism here is a recurring resonant phenomenon 
between the nonlinear wave propagation of the superfluid 
and the defect. The resulting scenario is plotted in Fig 1 
as a solid line, in contrast to the very general notion of 
the absence of superfluidity for v > Vc (dashed line). 

Superfluid motion in presence of defects has been ex¬ 
tensively investigated in ultracold atoms: superfluidity 
can be probed by stirring a laser beam [4, 10-12] and the 
critical velocity has been measured [4]. Experiments on 
superfluid motion have been performed also with mov¬ 
ing optical lattices [13], in toroidal geometries [14, 15], 
with ultracold fermions near unitarity [5] and in two- 
dimensional Bose systems [16]. 

To illustrate the existence of supercritical flow and 
the associated arising phenomena we choose defects of 
rectangular shape in the one-dimensional (Id) flow of a 
superfluid, whose behaviour is governed by the Gross- 
Pitaevskii equation. We consider a homogeneous system 
with stationary flow at velocity v in an initially flat po¬ 
tential, in which a rectangular defect is then introduced 
to study the dynamical response. From a theoretical per¬ 
spective, a challenging point is to define the transmission 
and reflection coefficients since the usual definitions used 
for linear matter waves [17, 18] do not apply anymore; 
in fact the superposition principle of an incoming and 
a reflected wave is no longer valid, and bound states in 
the defect can be present due to the interaction term 
[19]. For the interacting gas one can quantify the trans¬ 


mitted part of an incident wavepacket [20-26] or char¬ 
acterize the breakdown of superfluidity caused by the 
drag exerted by a matter wave on an obstacle [27]. In 
presence of a 5-like potential moving at supersonic veloc¬ 
ity the stationary wave patterns have been studied [28]. 
When the barrier is rectangular, solutions of the time- 
independent Gross-Pitaevskii equation can be written in 
terms of Jacobian functions [29, 30]: for a rectangular 
well the current-phase relation for subsonic motion can 
be determined [31, 32] and when neglecting the mean- 
field interaction outside the potential well, it is possible 
to analytically calculate the transport properties of the 
system in terms of incoming and outgoing waves and reso¬ 
nances and bound states are obtained in closed form [33]. 
Understanding matterwave propagation in the presence 
of tailored defect potentials is important for the study of 
analogue gravity models and acoustic Hawking radiation 
in Bose-Einstein condensates [34]. 

It is well known that in the non-interacting regime 
the transmission coefficient across a square potential, as 
found by solving the time-independent linear Schrodinger 
equation, reaches exactly unity for specific values of 
the momentum of an incident plane wave (Ramsauer- 
Townsend resonances). In this work we quantify the cre¬ 
ation of excitations, given the geometric parameters of 
the defect, in relation to the flow velocity v and the in¬ 
teraction strength g in the gas. We find that the resonant 
behaviour has a counterpart in the nonlinear regime and 
superfluid motion across a defect is possible in spite of the 
Landau criterion predicting growth of excitations. The 
velocities at which this excitation-free supercritical flow 
occurs are shifted and continuously connect to the veloc¬ 
ities of unit transmission in the linear case as the mean- 
field interaction parameter is varied. This can be viewed 
as generalised Ramsauer-Townsend resonances that con¬ 
tinuously connect the non-interacting case to the attrac¬ 
tively and repulsively interacting regimes. 

Given the explicit arguments of the Landau criterion 
it is important to clarify what is meant by supercriti¬ 
cal superfluid flow: i) K supercritical solution of the dy¬ 
namical equations exists for specific values of the super¬ 
fluid’s momentum and system parameters (i.e., interac¬ 
tion strength and geometric properties of the defect) and 
the flow of such a solution is stable under small pertur¬ 
bations. a) The perturbations of the propagating flow 
created by the defect is bounded in time (at least for 
experimental timescales) and no new excitations are pro¬ 
duced when the barrier is completely ramped on. Hi) The 
superfluid flow is physically inducible and accessible. We 
have found that all three conditions are fulfilled for the 
supercritical flow discussed here. 

We have verified that the production of excitations 
at the supercritical excitation-free points is bounded 
in time, that the flow is stable under small perturba¬ 
tions and for times larger than the typical experimental 
timescales of ultracold atom experiments with Id Bose 
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gases, and that the obtained findings do not crucially 
depend on the ramping time of the well/barrier. 

The Gross-Pitaevskii equation we consider is 

where is the condensate wavefunction, g is the 

one-dimensional nonlinear coefficient [35], and the po¬ 
tential is chosen to be 

, ( 2 ) 

[ 0 otherwise 

where d is the width and Vq the strength of the de¬ 
fect, which is ramped on at time t = 0 with a speed 
parametrized by a. At time t = ^barrier = 1.5a the value 
of the defect is 0.98Vb, so it is practically almost com¬ 
pletely turned on. The initial state (when the barrier is 
absent) is a plane wave with momentum k and velocity 
V = hk/m: = 0) = and density n = jV’oP- 

Since for t = 0 the defect is absent, il^{x,t = 0) is a so¬ 
lution with momentum k of the time-independent Gross- 
Pitaevskii equation. The idea of ramping on the defect is 
to adiabatically lead the system towards possible superc- 
tical superfluid solutions for specific values of the mo¬ 
mentum k. 

We consider periodic boundary conditions on a domain 
with finite length 2L d (with x G [—L,L]), and we 
check that for the considered times the effects of raising 
the barrier do not propagate to the boundaries; thus, the 
results are independent of the specific choice of bound¬ 
ary conditions. For the numerical solution we used the 
projected fourth-order Runge-Kutta method in the inter¬ 
action picture [36] , which we confirmed to be very stable 
for all times within the considered range. 

To quantify the production and growth of excitations 
we found convenient to introduce the time-dependent 
quantity D{t) (which we refer to as the disturbance) de¬ 
fined as 

= y dx{\i;{x,t)\'^-\'ilj{x,t = 0)\'^Y , (3) 

where the integral is calculated over the region where the 
initial wave propagation is directed towards the defect. 
j'0(a;,t = 0)]^ is the density n. Of course, if the defect is 
absent D = 0 . 

Our results are summarized in Figs. 2-4. Fig. 2 shows 
that essentially no excitations are produced for velocities 
below the sound velocity c = ^/gn/m (here Vc ~ vl = c). 
As expected, excitations are produced for v > Vc, as time 
progresses D{t) increases, and for large velocities D{t) 
becomes smaller. However, there are velocities v > Vc, 
for which the production of excitations is inhibited. As 
shown in the inset of Fig. 2, close to these points the 
production of excitations is bounded in time. Panels (a). 



FIG. 2. Disturbance 2) vs u/uc at various equidistant times 
(growing from the bottom of the figure), where Uc is equal to 
the sound velocity c. Darker lines indicate longer times. The 
line highlighted in red corresponds to t = tbarrier. At specific 
values of v, the disturbance does not grow with time after 
the barrier has finished rising. The insets illustrate this phe¬ 
nomenon in detail for the first minimum. Panels (a),(b),(c) 
show the density n{x) = lipix, to)|^ (in blue) and the potential 
V{x,to) (in green) at time to = 2.2tbarrier for three different 
initial values of v corresponding to the minimum and max¬ 
imum points for D indicated as “a”, “b”, “c” in the main 
figure. The initial condition is a plane wave travelling with 
momentum hk = mv towards the right. The simulation pa¬ 
rameters are g = 15 x (ft^/2md), a = 3 X (2md^ /fi) and 

Fo = 2 X /2md^Y 


(b), (c) of Fig. 2 show the density at a fixed time larger 
than tbarrier. Away from the minima phonons are emit¬ 
ted (Fig. 2b) while at the minima a stationary breathing 
state forms inside the defect (Fig. 2a and Fig. 2c) [27]. 
We checked that these results neither depend on the par¬ 
ticular choice of the measure of disturbance considering 
other quantifiers of the excitations, nor on the periodic 
boundary conditions or the choice of L. These findings 
also do not critically depend on the value a. This is true 
for typical experimental values of a as well as for a —>■ 0 , 
where small fluctuations of S)(t) occur for specihc choices 
of quantifiers of the disturbance. We finally observe that 
our results hold for a wide range of interactions, includ¬ 
ing large values of 5 > 0 , so that the validity spans the 
regime where the healing length is smaller than the defect 
width to the regime where it is larger. For small negative 
values of g the resonances are still present, but when —g 
becomes large then collapse and instabilities are observed 
as expected. 

As shown in Fig. 3, there is a clear difference between 
the growth of 23 (t) at the resonant and non-resonant ve¬ 
locities. At resonant velocities, excitations are produced 
exclusively during the ramping of the defect and the dis¬ 
turbance is afterwards constant. For non-resonant veloci¬ 
ties the disturbance grows linearly with time, which leads 
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Time {2m<f/h) 


FIG. 3. Disturbance T) as a function of time for the initial 
velocities v corresponding to the minima and maximum indi¬ 
cated in Fig. 2 (top) - the bottom plot shows the correspond¬ 
ing time dependence of the barrier strength. The background 
shading indicates the rate of change of the barrier strength. 
For times t > tbarrier, where the barrier strength has settled, 
D has a linear dependence on time. The creation of excita¬ 
tions can be characterised by the slope of the disturbance D 
at large times and it is extremely small at the minima (see 
central right panel in Fig. 4 for details). 



FIG. 4. Resonant momenta (with murea = lifcres). Centre- 
right: disturbance growth rate J) as a function of kd for four 
different values of interaction: gnjVo =-0.5, 0, 0.5, 7.5 (blue, 
green, red, purple lines, respectively). Left: Resonant mo¬ 
menta fcresd as a function of gn/Vb for the first four excitation- 
free points. The solid red line is the perturbative prediction 
(Eq. 4), the blue points indicate numerically calculated val¬ 
ues and the green line shows the sound velocity c = -y/ gn/m. 
Bottom-right and top-right: detailed view of the first and 
fourth resonances, respectively. 


us to quantify the growth of excitations by computing the 
time derivative of 23 (t) for t > ^barrier (and checking that 
the obtained value does not depend on the computation 
interval). The analysis shows that 2) is extremely small 
at the resonant velocities (being suppressed by at least 4 
orders of magnitude with respect to non-resonant veloc¬ 
ities). 2) is also very small for u < Uc in agreement with 
the Landau criterion. In the inset of Fig. 4 we plot the 
values 23 as a function of the velocity for four interaction 
strengths. By performing the same analysis we can re¬ 
construct the behaviour of the resonant velocities Ures as a 
function of g: the results are plotted in Fig. 4 for a barrier 
(positive Vb)- With respect to the non-interacting limit 
g = 0 the shift is positive (negative) for repulsive (attrac¬ 
tive) interactions g > 0 {g < 0). Similar results hold for 
a potential well (Vb < 0), then with a negative (positive) 
shift for repulsive (attractive) interactions. Our results 
are plotted in Fig. 4 together with a multiple-scale an¬ 
alytical derivation valid for small g, predicting values of 
^barrier (defined below) corresponding to a total transmis¬ 
sion across the barrier [37]. These momenta are given by 

^barrierC? = UTT + 6, where /Cbarrier = 

velocity v = hk/m) and 

- 3a(kd + mr) 

8 <“> 


with a = 2mgn/v^. Our result differs by the factor (—1)" 
from Eq. 24 of [37] . The analytical predictions for small 
g match the numerical data well for the higher excitation- 
free points, but less so for the lower resonant velocities. 
The ratios of the slopes of Uj-es as a function of g between 
the analytical and the numerical results are 0.47, 0.83, 
0.92 and 0.95 for the first four excitation-free points in 
increasing velocity order. 

An experimental setup to test these results can be im¬ 
plemented with ultracold Bose gases trapped on an atom 
chip. A Id quasi condensate can fill a potential tube 
of several hundred micron length [38] with a very small 
potential variation along the tube and correspondingly 
homogeneous density. The gas can be set in motion at 
a controlled velocity by removing the residual longitudi¬ 
nal confinement and applying a short pulse of a magnetic 
gradient in the same direction. Velocities on the order of 
and exceeding a typical sound velocity of c ~ 1 mm/s 
can be achieved straightforwardly. By applying currents 
to microwires on the chip, a magnetic defect can be pro¬ 
duced and controlled. Its geometric shape can be tailored 
with a resolution given by the atom-surface distance zq. 
Here it is critical that Zq is on the order of single microns 
in order to distinguish the individual excitationless res¬ 
onances from the intermittent regimes of fast excitation 
growth. The excitation behaviour can be probed by vary- 
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ing the velocity v of the gas’s motion or by varying the 
final amplitude of the defect Vq at fixed v. The difference 
in Vo for the first two excitation-free points is expected to 
be AVb « 1^32 = h X 6.9 kHzx/xm^Jj- for the example 
of ®^Rb, so that d > zq needs to be sufficiently small to 
maintain AVq « /i « lixl kHz, where /x is the chemical 
potential of the repulsively interacting gas. Appropriate 
rise times of the defect are on the order of milliseconds 

(^barrier > Ims) for Zq > 

In conclusion, we have studied the propagation of mat¬ 
ter waves across defects of rectangular shape starting 
from a stationary flow solution in the defect-free case 
with velocity v and ramping on the defect: for velocities 
smaller than a critical velocity there is no production 
of excitations (with Vc very well approximated by the 
Landau critical velocity vl in our one-dimensional case). 
For a set of supercritical velocities v > Vc the growth of 
excitations is fully suppressed, contrary to the generic ex¬ 
pectation as a consequence of the Landau criterion. For 
these velocities, we find the production of excitations to 
be bounded in time and to stop entirely when the defect 
is completely turned on. The flow is stable for times far 
in excess of any typical experimental timescale in the do¬ 
main of ultracold atoms and elongated Bose condensates 
in particular. The obtained findings do not crucially de¬ 
pend on the ramping time of the defect. Such excitation- 
free supercritical velocities are present both for wells and 
barriers, and for repulsive and (small) attractive interac¬ 
tions. We observe that even though in the nonlinear case 
bound states and bifurcation effects are expected, our 
protocol allows us to access the excitation-free points in 
a clean way, not depending on the ramping time. 

The obtained excitation-free supercritical velocities are 
the nonlinear counterpart of the velocities having unit 
transmission in the linear Schrodinger case (Ramsauer- 
Townsend resonances) and are due to the resonance be¬ 
tween the length scale associated with the matter wave 
momentum (^ 27r/fc) and the length scale of the defect. 
The shift from the resonance is positive (negative) for 
repulsive (attractive) interaction in the case of barrier 
defects, and vice versa for well defects. Beyond this sur¬ 
prising proof of the presence of such excitation-free super¬ 
critical velocities for the paradigmatic case of rectangular 
defects, we expect that they exist in general whenever 
a defect can be characterised by a well defined length 
scale, e.g. for trapezoidal defects or two delta-peaked po¬ 
tentials. The steeper the defect is at its edges, the more 
robust will the inhibition of excitations be in the vicinity 
of a set of supercritical velocities. 

Motivated by our work it will be interesting to study 
general criteria of existence and stability of supercritical 
solutions of the nonlinear Schrodinger equation in pres¬ 
ence of defects with a well defined length scale. On the 
other hand, intriguing possibilities arise from utilising 
supercritical points in measurement devices based on 
superfluids and superconductors. Tuning a barrier to 


a supercritical flow resonance would facilitate precise 
determination of unknown external parameters affecting 
the flow velocity such as rotation and performing 
selective measurements at supercritical velocities. 

Acknowledgements: Discussions with S. Fagnocchi, G. 
Dell’Antonio, P. Kevrekidis , G. Mussardo and A. Recati 
are gratefully acknowledged. A.T. acknowledges support 
from the Italian PRIN “Fenomeni quantistici collettivi: 
dai sistemi fortemente correlati ai simulator! quantistici” 
(PRIN 20I0_20I0LLKJBX). A.P.M. acknowledges sup¬ 
port from CONACyT. Support from the EU-FET Proac¬ 
tive Action (Grant 601180 “MatterWave”) and EPSRC 
(Grant EP/I017828/1) is also acknowledged. 


* Peter.Kruger@nottingham.ac.uk 

[1] M. Tinkham, Introduction to superconductivity (New 
York, McGraw-Hill, 1996). 

[2] A. J. Leggett, Quantum liquids: Bose condensation and 
Cooper pairing in condensed-matter systems (Oxford, 
Oxford University Press, 2006). 

[3] D. R. Allum , P. V. E. McClintock, and A. Phillips, Phil. 
Trans. R. Soc. A 284, 179 (1977). 

[4] C. Raman, M. Kohl, R. Onofrio, D. S. Durfee, C. E. 
Kuklewicz, Z. Hadzibabic, and W. Ketterle, Phys. Rev. 
Lett. 83, 2502 (1999). 

[5] D. E. Miller, J. K. Chin, C. A. Stan, Y. Liu, W. Setiawan, 
C. Banner, and W. Ketterle, Phys. Rev. Lett. 99, 070402 
(2007). 

[6] E.M. Lifshitz and L.P. Pitaevskii, Course of theoretical 
physics, V.9: Statistical physics, (London, Pergamon, 
1980). 

[7] G. Baym and C. J. Pethick, Phys. Rev. A 86, 023602 

( 2012 ) 

[8] S. Finazzi, F. Piazza, M. Abad, A. Smerzi, and A. Recati, 
arXiv:1409.8068 

[91 S. laneselli, C. Menotti, and A. Smerzi, J. Phys. B 39, 
S135 (2006). 

[10] R. Onofrio, C. Raman, J.M. Vogels, J.R. Abo-Shaeer, 

A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett. 85, 
2228 (2000). 

[11] P. Engels and C. Atherton, Phys. Rev. Lett. 99, 160405 
(2007). 

[12] T.W. Neely, E.G. Samson, A.S. Bradley, M.J. Davis, and 

B. P. Anderson, Phys. Rev. Lett. 104 , 160401 (2010). 

[13] J. Mun, P. Medley, G.K. Campbell, L.G. Marcassa, D.E. 
Pritchard, and W. Ketterle, Phys. Rev. Lett. 99, 150604 
(2007). 

[14] A. Ramanathan, K.C. Wright, S.R. Muniz, M. Zelan, 
W.T. Hill, C.J. Lobb, K. Helmerson, W. D. Phillips, and 
G. K. Campbell, Phys. Rev. Lett. 106 , 130401 (2011). 

[15] S. Moulder, S. Beattie, R.P. Smith, N. Tammuz, and Z. 
Hadzibabic, Phys. Rev. A 86, 013629 (2012). 

[16] R. Desbuquois, L. Chomaz, T. Yefsah, J. Lonard, J. 
Beugnon, C. Weitenberg, and J. Dalibard, Nature Phys. 
8, 645 (2012). 

[17] C. Cohen-Tannoudji, B. Din, and F. Laloe, Quantum me¬ 
chanics (New York, Wiley, 1977). 



6 


[18] R. Gilmore, Elementary quantum mechanics in one di¬ 
mension (Baltimore, Johns Hopkins University Press, 
2004). 

[19] P. Leboeuf and N. Pavloff, Phys. Rev. A 64, 033602 

( 2001 ). 

[20] J. S. Stiessberger and W. Zwereer, Phys. Rev. 62, 
061601(R) (2000). 

[21] T. Paul, K. Richter, and P. Schlagheck, Phys. Rev. Lett. 
94 , 020404 (2005). 

[22] T. Paul, M. Hartung, K. Richter, and P. Schlagheck, 
Phys. Rev. A 76 , 063605 (2007). 

[23] A. M. Kamchatnov and L. P. Pitaevskii, Phys. Rev. Lett. 
100 , 160402 (2008). 

[24] A. M. Leszczyszyn, G. A. El, Yu. G. Gladush, and A.M. 
Kamchatnov, Phys. Rev. A 79 , 063608 (2009). 

[25] G. Watanabe, F. Dalfovo, F. Piazza, L. P. Pitaevskii, and 
S. Stringari, Phys. Rev. A 80 , 053602 (2009). 

[26] F. Piazza, L. A. Collins, and A. Smerzi, New J. Phys. 
13 , 043008 (2011); J. Phys. B 46 , 095302 (2013). 

[27] N. Pavloff, Phys. Rev. A 66, 013610 (2002). 

[28] T.-L. Horng, S.-C. Gou, T.-G. Lin, G. A. El, A. P. Itin, 


and A. M. Kamchatnov, Phys. Rev. A 79 , 053619 (2009). 

[29] B. T. Seaman, L. D. Carr, and M. J. Holland, Phys. Rev. 
A 71, 033609 (2005). 

[30] K. Rapedius and H.J. Korsch, Phys. Rev. A 77 , 063610 
(2008). 

[31] A. Baratoff, J. A. Blackburn, and B. B. Schwartz, Phys. 
Rev. Lett. 25, 1096 (1970). 

[32] F. Piazza, L. A. Gollins, and A. Smerzi, Phys. Rev. A 
81, 033613 (2010). 

[33] K. Rapedius, D. Witthaut, and H.J. Korsch, Phys. Rev. 
A 73 , 033608 (2006). 

[34] See an updated review in C. Barcelo, S. Liberati, and M. 
Visser, Living Rev. Rel. 8, 12 (2005). 

[35] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[36] A. Norrie, Ph. D. Thesis (University of Otago, 2005) 

[37] H. A. Ishkhanyan and V. P. Krainov, Phys. Rev. A 80, 
045601 (2009). 

[38] P. Kruger, S. Hofferberth, 1. E. Mazets, 1. Lesanovsky, 
and J. Schmiedmayer, Phys. Rev. Lett. 105, 265302 
( 2010 ). 



